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Abstract. Fireballs with high energy per baryon rest mass (tj > Tj* ~ 400) contain a relatively 
slow neutron component. We show here that in this situation the thermal history of fireballs is very 
different from the standard adiabatic cooling. 



INTRODUCTION 

Baryonic matter ejected in a GRB explosion is partially composed of free neutrons (see 
i\\ and [2] for detailed calculations and references). This has a strong impact on the 
external blast wave and changes the afterglow mechanism at radii r ^ 10^^ cm [3]. We 
here investigate the dynamical effect of neutrons on the early evolution of a relativistic 
fireball, prior to the afterglow phase. 



FIREBALL ACCELERATION AND NEUTRON-PROTON 

DECOUPLING 

At the initial stage of the fireball acceleration, the neutrons are coUisionally coupled 
to the protons and accelerated by radiation pressure at the same rate. If the energy-per- 
baryon, 7], is high enough, neutrons decouple before the acceleration stage is completed. 
This happens if 
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where ^ is the ratio of neutron and proton densities, Rq = 10^ cmi?07 is the initial size of 
the fireball, and L = 10^^ergL52 is the equivalent isotropic luminosity of the explosion. 
Beyond the decoupling radius. 
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neutrons coast with a constant Lorentz factor 



^n^^^^nJ^X^'. (3) 
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The optically thick proton fireball continues to accelerate as Fp ~ r/i?o until the 
internal energy equals the proton rest-mass energy and acceleration is no longer efficient. 
If the fireball remains optically thick during the whole acceleration stage, its Lorentz 
factor saturates at 

( 1 v<'n*, 

and the acceleration stage ends at radius 

Rs^rpRo. (5) 



ADIABATIC COOLING 

We first consider the thermal history of the fireball without neutrons. Electrons, protons, 
and radiation maintain a common temperature T in the early dense fireball via Coulomb 
collisions and Compton scattering, and T decreases adiabatically to very low values. 
During such adiabatic expansion, the photon-to-baryon ratio (j) = ny/nt, — const ~ 10^, 
and radiation completely dominates the internal energy and pressure. Therefore, T 
decreases according to the radiation adiabatic law with index 7= |. This continues until 
the fireball becomes transparent to radiation at the photosphere radius R^. Assuming 
R^ > Rs, we have 

R 

here Oj is Thomson cross section. After becoming transparent the plasma is still tracking 
the temperature of the (freely streaming) photons, which is constant in the plasma frame 
if the fireball coasts with a constant Lorentz factor. The electrons decouple thermally 
from radiation only when the Compton timescale {tc = Tp{3mec/AUrad<^T) where Urad 
is the radiation energy density) exceeds the expansion timescale R/c. This happens at 
radius 



Rer= V-773 • 

The protons decouple from the electrons when the Coulomb timescale tgp ~ 
17 T^/^Hg^Fp s exceeds R/c. The corresponding radius of e-p decoupling is 

RepC.{6S7:mpcY^,[^), (8) 



where Ts = T^/Tp is the temperature at the saturation radius, and Tq is the initial 
temperature of the fireball. The protons decouple thermally from the electrons before 
the electrons decouple from radiation if Rgp < Rej. One can show that Rgp/Rey ~ 

4 X IQ-^T^^'^{R^/Rsfl^ (with T, = Tsi x 10^ K). At radii R > mm{Rey,Rep), the 
protons are decoupled from radiation and cool adiabatically with index 7 = f ■ 

Thus, in the absence of neutrons, the fireball cools down adiabatically as T <x 
during the acceleration stage and as T oc r^'^/^ during the subsequent coasting stage 
(Tp ~ r]) up to the transparency radius. After thermal decoupling from radiation the 
protons cool as Tp r^^/^. The thermal history of a pure proton fireball is shown by 
dot-dashed curve in Fig.[T] For most of the early evolution, the fireball is a cold coasting 
outflow. The presence of neutrons and their /3 -decay change this picture. 



HEATING BY PROTON-NEUTRON COLLISIONS 

Protons in the fireball are continuously heated via proton-neutron collisions The 
coUisional heating reaches its peak atR^ Rnp'. aiR< Rnp the n-p collisions are frequent 
but the relative velocity of the neutron and proton components is small (/3„p ~ < 1) 

and ai R> Rnp the relative velocity is relativistic /3„p ~ 1 however the collisions are 
rare. The coUisional heating, thus, peaks at Rnp where the relative velocity becomes 
comparable to the speed of light and the rate of collisions still allows an efficient transfer 
of energy on a dynamical timescale. There are two sinks of heat gained by protons via 
n-p collisions: adiabatic cooling and Coulomb scattering off electrons. The competition 
between heating and cooling shapes the first peak in the proton temperature profile 
(Fig. 1). At i? = Rep the energy transfer from the protons to the electrons becomes 
inefficient on the expansion timescale R/c. Then most of the heat gain by protons 
remains stored in the proton component, not given to the electrons and radiation. 



FIREBALL HEATING AND DECELERATION BY DECAYED 

NEUTRONS 

In a fireball with 7] > 7]*, the neutrons have a smaller momentum than the protons after 
Rnp. Their decay products e and p exchange momentum with the ion fireball and heat 
it up. The heating of an ion medium by /3 -decay of neutrons moving with respect to the 
medium was discussed in [3]. The decay particles exchange momentum with the medium 
through the two-stream instability or because of gyration in a transverse magnetic field 
frozen into the medium. Since the neutrons decay gradually at all radii, this momentum 
exchange and heating take place continuously from the moment of decoupling. 
The heating rate of the fireball due to /3 -decay is 

dEh _ dMnpC^ 

^-(r../-!)^;-, (9) 



where F^g/ is the neutron Lorentz factor in the fireball rest frame, 



'-^ = 1^, (10) 



is the decayed neutron mass per unit radius, and 

^^3:^ 0.8x1016 cm (J^j (11) 

is the mean radius of decay. The heating by /3 -decay begins at R^p and significantly 
changes the thermal history of the fireball. It shapes the second peak in proton tempera- 
ture (see Fig.U]). 

As long as proton thermal energy is much below nipC^, adiabatic cooling is small 
compared to heating and the proton temperature increases linearly with radius, 

TpOcr. 

At a radius R^us the dissipated bulk kinetic energy equals the total energy of the fireball 
and Tp begins to decrease. Between R^is and Rp the /3 -decay heating balances the 
adiabatic cooling and 
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When the fireball reaches Rp the neutron component is exponentially depleted and 
the dissipation process switches off. The adiabatic cooling then leads to a power-law 
decrease in temperature. 



TpOcr-'/\ 



DISCUSSION 

We have shown here that the presence of a neutron component greatly affects the early 
dynamics of GRB fireballs with r] > r]* ~ 400. Contrary to the pure proton model, 
the fireball with a neutron component heats up significantly at R ^Rnp- In a popular 
GRB scenario, the bulk energy is partially converted into the observed 7— rays through 
shock waves inside the fireball at r ^ 10 cm lH. The /3 -decay process described here 
should affect the internal shocks. It decelerates the fastest portions of the inhomogeneous 
fireball and reduces the contrast of Lorentz factors, which should significantly reduce 
the dissipation efficiency of internal shocks. A detailed study of these effects is in 
preparation. 
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FIGURE 1. Proton temperature as a function of radius in a fireball with neutron-to-proton ratio ^ = 1 
(solid curve). The relevant radii R^p, R„p, Rr, Rs, and Rf^ are shown by vertical dashed lines (note that 
Rs = Ri in this case). For comparison, the dot-dashed curve shows the evolution of a pure proton fireball 
((^ = 0) with the same luminosity L = 10^^ erg s^', initial radius Rq = 10^ cm and final Lorentz factor 
Tp ~ 1040. 
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FIGURE 2. The bulk Lorentz factor of protons (solid fine) and neutrons (dot-dashed line) as a function 
of radius for the same parameters as in Fig.Q] 
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